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Abstract 

In the model where two massive scalar particles interact by the ladder exchanges of 
massless scalar particles (Wick-Cutkosky model), we study in light-front dynamics 
the contributions of different Fock sectors (with increasing number of exchanged 
particles) to full normalization integral and electromagnetic form factor. It turns 
out that two-body sector always dominates. At small coupling constant a <C 1, its 
contribution is close to 100%. It decreases with increase of a. For maximal value a = 
2tt, corresponding to the zero bound state mass, two-body sector contributes to the 
normalization integral 64%, whereas the three-body contribution is 26% and the sum 
of all higher contributions from four- to infinite-body sectors is 10%. Contributions 
to the form factor from different Fock sectors fall off faster for asymptotically large 
Q 2 , when the number of particles in the Fock sectors becomes larger. So, asymptotic 
behavior of the form factor is determined by the two-body Fock sector. 
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1 Introduction 



In field theory, the decomposition of the state vector in the basis of free fields 
with given momenta results in the concept of relativistic wave function in 
momentum space. The latter is the coefficient of the decomposition - the 
Fock component. The state vector is described by an infinite set of the Fock 
components, corresponding to different numbers of particles. 
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However, the models used in applications usually consider only few-body Fock 
components. The infinite set of the Fock components is truncated to one com- 
ponent only, with two or three quarks. The belief that a given Fock sector 
dominates is mainly based on intuitive expectations and "experimental evi- 
dences" rather than on field-theoretical analysis. Without this analysis a model 
remains to be phenomenological. Though one can perturbatively estimate the 
next Fock components (with one or even few extra particles), this says nothing 
about convergence of full Fock decomposition. 

The difficulty to take into account the many-body sectors is caused by the fact 
that it is equivalent to solving a true field-theory problem. Important part of 
this problem - finding the Bethe-Salpeter (BS) function [1] - is solved in the 
Wick-Cutkosky model [2]. In this case the BS equation is reduced to the one- 
dimensional one which can be easily solved numerically for arbitrary coupling 
constant. In weak and strong coupling limits there are explicit analytical solu- 
tions. This allows us to use Wick-Cutkosky model for non-perturbative study 
the Fock decomposition. 

In the Wick-Cutkosky model two spinless constituents with mass m interact 
by exchange by a massless scalar particle. Corresponding BS function is the 
"two-body" one in the sense that it depends on two four-dimensional particle 
coordinates. The term "two-body" here is rather slang than reflection of a 
real physical situation. It may be misleading. We emphasize that this does 
not mean that the "two-body" BS function describes a two-body system, i.e., 
a system with the state vector truncated to the two-body Fock component. 
On the contrary, the corresponding state vector contains infinite set of the 
Fock components with two massive constituents and 0, 1, 2, . . . , oo massless 
exchanged particles. 



Fig. 1. Feynman ladder graph with two exchanges. 



This can be explained as follows. In terms of the Feynman diagrams, the 
one-boson exchange kernel generates the ladder graphs only. One of these 
graphs, with two exchanges, is shown in Fig. 1. We will work in the light-front 
dynamics (LFD). In terms of the time-ordered graphs, in the light-front (LF) 
time, this Feynman diagram contributes to the two-, three- and four-body 
states (with two massive constituents and with zero, one and two exchanged 
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n=2 n=3 n=4 n=3 n=2 

Fig. 2. One of six time-ordered ladder graphs, generated by the Feynman graph 
Fig. 1. 

bosons). There are six such time-ordered graphs. One of these graphs is shown 
in Fig. 2. Since the BS equation deals with infinite series of ladder exchanges, 
it implicitly incorporates all the intermediate states with number of particles 
from two to infinity. The problem (solved in the present paper) is to "extract" 
from the BS solution the two-body contribution and total contribution of all 
other Fock components. Comparison of total contribution with the two-body 
one shows the contribution of the many-body sectors. 

Because of the ladder kernel, the Wick-Cutkosky model is still approximate. 
All the exchange particles, emitted by one and the same boson, must be and 
absorbed by other one. The graphs with particles emitted by different bosons 
in the same intermediate state are generated by the cross Feynman box which 
is absent in the Wick-Cutkosky model. Self-energy graphs are absent too. In 
spite of these restrictions, the Wick-Cutkosky solution does contain in the state 
vector an infinite sum of many-body Fock components. Namely this property 
is most important for our work. 

In Wick-Cutkosky model there are also some abnormal solutions which have 
no counterpart in the non-relativistic potential theory [2,3]. There was a dis- 
cussion [3], whether they are physical ones or not. This is a separate problem, 
not related to our study. We deal with the normal solution. In non-relativistic 
limit the two-body Fock component, which only survives in the state vector 
of normal solution, turns into the usual ground state wave function in the 
Coulomb potential. 

A truncated LF field theory was developed in [4]. Some estimations of the 
higher Fock sectors were done in [5]. Truncated LF Fock space decomposition 
was applied in [6] to nonperturbative study the large a QED. Recently, the 
Fock decomposition, incorporating first three Fock sectors, was used for non- 
perturbative renormalization in a scalar model [7] and, with two sectors, - in 
Yukawa model [8,9] and in gauge theory [9]. 
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In the present paper, we study in the Wick-Cutkosky model the contribution 
of the many-body sectors in electromagnetic form factor and compare it with 
the two-body contribution. Besides, we calculate the contribution to the nor- 
malization integral of the two- and three-body states. Subtracting them from 
1, we get contribution of all the states with n > 4. In this way we investi- 
gate the convergence of the Fock decomposition. Calculations are carried out 
nonperturbatively in full range of binding energy < B < 2m. 

Plan of the paper is the following. In Section 2 we define the Fock decomposi- 
tion (in the framework of LFD), and remind the definition of the BS function. 
In Section 4 the BS equation is solved. In addition to analytical solutions in 
limiting cases M — > 2m and M = 0, found in [2], we find, for small a, next 
correction of the order of a logo; and solve BS equation numerically for any a. 
In Section 5 the two-body contribution to the normalization integral and to 
the electromagnetic form factor is calculated for both small and large coupling 
constants. The form factor asymptotic at Q 2 — > is obtained. In Section 6 
the same is done for three-body contribution and in Section 7 for sum of all 
Fock sectors. Sect. 8 is devoted to numerical calculations for arbitrary coupling 
constant. Sect. 9 contains summary and discussion. Some technical details are 
given in Appendices A, B and C. 



2 Fock decomposition and BS function 



The state vector \p), eigenstate of full Hamiltonian, can be decomposed in 
terms of the eigenstates of free Hamiltonian. This results in natural formu- 
lation of the concept of a relativistic wave function in terms of the LF Fock 
expansion at fixed LF time, which is usually put to zero for a stationary state: 
uj-x = 0. The null four-vector uj (u> 2 = 0) determines the orientation of the LF 
plane; the freedom to choose uj provides an explicitly covariant formulation of 
LF quantization [10]. The Fock decomposition has the form: 



\p) = J2 / ip{h,...,k n ,p,ujT) 

n>2 J 



X 



[Y.kj-p-ur] 2(uj-p) dr J[ ^JW(fc.) |0> . (1) 



Here a) is the usual creation operator and = ym 2 + k 2 . All the four- 
momenta are on the corresponding mass shells: k 2 = m 2 , p 2 = M 2 , (ujt) 2 = 0. 
In QCD, the set of LF Fock state wave functions ip{ki, . . . , k n ,p,ur) repre- 
sents the ensemble of quark and gluon states possible when the hadron is 
intercepted at the light-front. In Wick-Cutkosky model this set represents 
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the states with two massive constituents +0, 1,2,... massless exchange parti- 
cles, forming two-, three-, four-body states, etc. The constituents are spinless 
and we consider the state \p) with total zero angular momentum. Therefore 
we omit the spin indices. The scalar variable r controls the off-shell con- 
tinuation of the wave function. From the point of view of kinematics, the 
four-momentum ojt can be considered on equal ground with the particle four- 
momenta fci, . . . , k n ,p. Being expressed through them (by squaring the equal- 
ity EI" h = p + ut), t reads: r = (M 2 - M 2 )/{2uj-p), where M 2 = (£™ ki) 2 . 

For system with zero total angular momentum the wave functions ifj are the 
scalar functions and they depend on a set of scalar products of the four- 
momenta k\, . . . ,k n ,p, ut with each other. One should choose a set of 3n — 4 
independent scalar products. A convenient way to choose these variables is 
the following. We define 

Xi = u-ki/u-p , Rt = hi- Xip (2) 



and represent the spatial part of the four-vector Ri = (Rio, Ri) as Ri = R^ + 

— * — * — * 

Rij_, where Rq is parallel to uj and Ri± is orthogonal to Q. In the standard 
version of LFD the difference J2iki — p is non-zero for the minus component 
only. In covariant formulation this means: J2i ki = p + ur. Therefore: J2i Ri± — 
0, J2i x i — 1- Since Ri-oj = RioU — Ri\yCu = by definition of Ri, it follows 
that R40 = \Ri\\\, and, hence, Rf ± = —Rf is Lorentz and rotation invariant. 

— * — * 

Similarly one can show that Ri±-Rj± = —RfRj. Therefore the scalar products 
Ri±-Rj± are also Lorentz and rotation invariants. In terms of these variables 
the invariant energy s = Mq = (J2i ki) 2 is given by 

s = J2 R ^ +m \ (3) 

Xi 



The variables Ri±, Xi are analogous to the well-known variables in the infinite 
momentum frame [11]. 

— * — * — * 

The wave functions are parametrizd as: ip = 4>(Ri±,Xi, R2±,x 2 ; ■ ■ ■ ; R n ±,x n ) 
and should depend on Xi, R 2 L and on Ri±-Rj±. In terms of these variables, the 
integral in Eq. (1) is transformed as 



d 3 L; 



i (2vr) 3 2e fci 

En I TT d 2 Ri± d.Xi 
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The state vector (1) enters the definition of the BS function: 

$(x 1 ,x 2 ,p) = (0\T((p(x 1 )<p(x 2 ))\p). (4) 



Since the function $(xi, x 2 ,p) depends on two four-dimensional coordinates 
x±,X2, one often refers to the function (4) as describing a two-body system. 
However, the state vector \p) is the full state vector. It contains all the Fock 
components, with the constituent numbers n = 2,3,4,.... The operators 
(p(x 1 ) : (p(x 2 ) in (4) are the Heisenberg operators (of the massive field). 

We need the BS function $(/c,p) in the momentum space. It is related to (4) 
as follows: 



$(x 1 ,x 2 ,p) = (2tt) 3/2 exp(-^-(a;i + x 2 )/2)$(x,p), 

&(k,p) = J <P(x,p) exp(ik-x) d 4 x, (5) 

where x = x± — x 2 and p = p\ + p 2 is the on-mass-shell momentum of the 
bound state (p 2 = M 2 ), p±,p 2 are the off-shell momenta {jp\ ^ p\ ^ m 2 ) and 
k = (px -p 2 )/2. 

Knowing the BS function (5), we can extract from it the two-body component 
[10]: 

i>fa,x) = tb^J 7 Hk+Mdp . ( 6 ) 

7r{u-p) J 



This relation is independent of any model. 



In the Wick-Cutkosky model, the BS function for the ground state 

with zero angular momentum L = has the following integral representation 
[2,3]: 

v ) = % _ 7 9m(z) dz 

1 ,P) x/4^_i (m 2 - M 2 /4 - k 2 - zp-k - ^e) 3 ' 1 } 



This representation is valid and exact for zero-mass exchange. The function 
gu{z) is determined by an integral equation [2,3] given below in Section 4. It is 
solved analytically in the limiting cases of asymptotically small and extremely 
large binding energy and numerically - for any binding energy. 
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Substituting (7) into (6), we find the two-body Fock component [12,13]: 



nR±,x) = — — -a. (8) 

2 V^[R± +m 2 - x{l - x)M 2 



When we express electromagnetic form factor through the BS function, we 
do not make any truncation of the Fock space. We get sum of contributions 
of all the Fock components. On the other hand, having found the two-body 
component (8), we can calculate its contribution to form factor and to the 
full normalization integral. In addition, below we will calculate explicitly the 
three-body contribution. Comparison of full contribution with the two-body 
ones allows to find the many-body sector contribution, for arbitrary values of 
coupling constant. Comparison of two-body, three-body and full contributions 
allows to see, how the Fock decomposition converges. In the present paper we 
will realize this program. 

We will use the BS function (7) and the LF two-body wave function (8) to cal- 
culate full form factor (and normalization integral) and two-body contribution 
to them. For these calculations we need to know the function gu( z )- 



3 Normalization condition for gu(z) 



Since the state vector \p) in definition (4) of the BS function is normalized, the 
BS function does not contain any arbitrary factors and, hence, is also prop- 
erly normalized. The normalization factor cannot be found, of course, from the 
homogeneous BS equation, but is determined by a special normalization con- 
dition [3]. It is equivalent to the condition based on the charge conservation [3], 
which means that Ff u u(0) = 1. This condition determines the normalization 
of g M (z). 




P k p' 

Fig. 3. Feynman diagram for the EM form factor. 
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The electromagnetic form factor Ff u u(Q 2 ) (Q 2 = —q 2 > and q = p' — p) 
is shown graphically in Fig. 3. The vertices at the left- and right-hand sides, 
being multiplied by the propagators, are the BS functions (5). Therefore the 
electromagnetic vertex is expressed in terms of the BS function as [10] 



{p + pyF full {Q 2 ) (9) 
/ (2^F (P + P'~ 2k ^ {m 2 -k 2 )^{^p-k,p^{^p' -k,p^. 



We multiply both sides of (9) by (p + p')^, insert (7) into (9), use the Feynman 
parametrization: 

1 } 30w 2 (l -u) 2 du 



° 3&3 / (au + bil-u))*' 

integrate (using Wick rotation) over four-vector k and replace the variables 
z — 2x — 1, z' — 2x' — 1. In this way, at Q = we find the normalization 
condition: 

F fM0) = -35^3 / 9m(2x -l)dx J g M {2x' - 1) dx' (10) 



} 2„ \2 A ((6e-5)m 2 + 2e(l-OM 2 ) 
x / u (1 u) 2 du ± — J - = 1, 

J ( m 2_ e(1 _ 0M 2) 

where £ = xu + x'{l — u). 

4 Solution of Wick-Cutkosky equation 

The integral equation for the function gu{z) for the ground state reads [2,3]: 



i 

a 



9m{z) = — I K(z,z')g M (z')dz' (11) 



-i 



with the kernel: 



K(z,z') = 



m 2 



m 2 - \{l - z' 2 )M 2 



(! ~ z ) a (y ^ , (j + f) n(J ^ 



.(12) 



Here a = g 2 /(167im 2 ) and g is the coupling constant in the interaction Hamil- 
tonian H mt = —g(p 2 (x)x{x). In nonrelativistic limit the interaction is reduced 
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to the Coulomb potential V(r) = —f. The function 9m{ z ) is defined in the 
interval — 1 < z < 1 and it is even: 9m(— z ) = 9m( z )- 

At small binding energy B = 2m — M, and, on the contrary, at extremely 
large binding energy equal to 2m, the solutions 9m{. z ) are found explicitly. At 
asymptotically small binding energy (corresponding to a — > 0) it has the form 
[2]: 

9M^2m(z) = N M ^ 2m (l - \z\) (13) 

and the binding energy is the non-relativistic one in the Coulomb poten- 

2 

tial: B = m f-. The normalization constant Nu^im is found by substituting 
9M^2m( z ) into the normaization condition (10). Calculation gives: 

N M ^ 2m = 8V2na 5/2 m 3 . (14) 



The higher Fock sectors are generated by extra exchanges which contain extra 
degrees of a. They are omitted in the solution (13). Therefore, to analyze the 
many-body contributions, we should take into account next ~ a correction to 
Eq. (13). To find it, we use the method of paper [14]. The correction contains 
both the terms ~ a and ~ a log a. We keep the leading ~ a log a only. The 
details are given in Appendix A. The solution of Eq. (11) obtains the form: 



9m{z) 



N 



a 



I - |.:| f — (1 + \z\)\og(z 2 + a 2 /A) 

Z7T 



(15) 



Corresponding binding energy B = 2m — M is given by [14]: 



B 



ma 2 ma 3 



log-. 

7T a 



(16) 



The normalization factor iV in (15) (still found from the condition (10)) now 
reads: 



N = 8V27ra 5/2 m 3 



5a _ 
1 H log a 



(17) 



In the opposite limiting case B = 2m, M = 0, which is achieved at a = 2n, 
the solution has the form [2,3]: 



9m=o(z) = N M=0 (1 - z 2 ) . 



(18) 



The normalization coefficient reads: 
N M =o = QV307r 3/2 m 3 . 



(19) 
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The latter case is especially simply checked. We substitute Eq. (18) with the 
normalization (19) into Eq. (10) and find: 

l 

F,,„(0) = 540/^(1 - uMl - - x<)(5 - 3* - 3*<) d U d*d*< = 1. 



This justifies the value (19). 

The BS functions corresponding to the small and large binding energy solu- 
tions (13) and (18) are given explicitly in Appendix B. 

For arbitrary bound state mass M we solve equation (11) numerically and 
find corresponding function gu[z) and coupling constant a. The interval 
— 1 < z < 1 is split in 50 equal subintervals. In every subinterval the function 
gu^z) is represented by a sum of the quadratic spline functions. We substi- 
tute a given initial <?m(z) in the right-hand side of (11), represent the result 
of numerical integration again through the spline functions and iterate. For 
M — * 2m and M — > we numerically reproduce the analytical solutions (13) 
and (18) correspondingly. The iterations converge very quickly. For example, 
for M = 0, even if we start iterations with the function (13), corresponding to 
the opposite limit M — > 2m, after a few iterations we reproduce the true solu- 
tion (18). The accuracy of seven digits (for example) in the function gM=o(z) 
is achieved after 8 iterations. 

In Fig. 4 the dependence of M on a is shown. In all the figures we use the 
units of to, i.e. we put to = 1. One can see that M turns into zero at a ~ 6.28, 
that agrees with the analytical solution a(M = 0) = 2n. 




1 2 3 4 5 6 7 

a 



Fig. 4. Bound state mass M vs. coupling constant a. In all the figures we put m = 1. 

To compare the form of the functions gu(z) vs. z for different M's, we, at first, 
normalize all Qm{z) 's by the condition <?m(0) = 1 and show them in Fig. 5 for 
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1 


9m( z > 


M=0 


/ /// 08 ~ 




\ M=1.9 


/ 06 ~ 




\ \ A M — 1 QQ 

\ \/ \ ivi — i .yy 
\\X\ M=1.9999 


/ /// 04 " 






/// 02 ~ 






, f , 1 , 0.0 


1 


i 1 l 1 



-1.0 -0.5 0.0 0.5 1.0 



Fig. 5. Functions gu{z) normalized by #m(0) = 1. 

the values M = 0, 1.9, 1.99, 1.9999. It is seen that the form of Qm^z) does not 
depend significantly on M. For very different values M = and M = 1.9 the 
curves of Qm{z) are rather close to each other. All the curves for < M < 1.9 
are between these two curves. We emphasize that with the function gu{z) 
normalized by <?m(0) = 1 the form factor Ff u u(0) is not normalized to 1. 

The functions Qm{z) normalized by the condition Ff u u(0) = 1, Eq. (10), are 
shown in Fig. 6 for the values M = 0, 1, 1.5, 1.75, 1.9, 1.9999. They strongly 
depend on M. This dependence comes mainly from the normalization factor. 
In the scale of Fig. 6 the function gM( z ) at M = 1.9999 is indistinguishable 
from the z-axis. 

Now we are in position to calculate form factors and two- and three-body 
contributions to the normalization integral. 



5 Two-body contribution 

5. 1 Form factor 

The two-body contribution F 2 b{Q 2 ) to form factor Ff u u(Q 2 ) is shown graphi- 
cally in Fig. 7. This diagram corresponds to the time-ordered (in the LF time) 
graph technique. The latter exists in a few versions [11,10,15], giving equiva- 
lent results. Dealing with the scalar particles, we will use the Weinberg rules 
[11] (given in Appendix C). Though the expression for two-body form factor in 
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Fig. 6. Functions gu(z) normalized by Ff u u(0) = 1. 




Fig. 7. LF diagram for the two-body contribution to the EM form factor. 

terms of an overlap of the LF wave functions is well-known [10,15,16,17,18,19], 
we derive it in Appendix C. Form factor reads: 



1 f -> -> -> d R 

F 2 b(Q 2 ) = -j^-p J ^(-Ri±, a:i)^(-Ri± ~ Xl ^^ x ^2x^ 



d 2 i?i_L dxi 



(20) 



where Q\ = Q 2 . Substituting in (20) the wave function ip(R±,x) determined 
by Eq. (8) and combining two i?u_-dependent denominators by means of the 
formula 

1 f 6u(l — u) du 



[ 

° 2&2 { (au + b(l-u)) 



4 - 
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we can easily integrate over R±±: 

1 1 

2> t : ' 



F 2b (Q 2 ) = 1 / / x ( 1 ~ x ) 9m(2x - 1) dx u(l - u) du ^ 



{ (u{l - u)x 2 Q 2 + m 2 -x{l-x)M 2 ) 
Calculating then the integral over u, we finally obtain: 



F 2b {Q 2 ) 



2% 3 g e 



o 



4 7 (1 + 7 ) log + (1 - 27)^1 + 47 

(l-x)g 2 M (2x-l)dx 

x 5 7 ( 1 + 4 7 )5/2 ' ^ > 



where 



7 = 



m 2 — x(l — x)M 



2 



x 2 Q 2 

5.2 Normalization integral 

5.2.1 Small binding energy (M — > 2m) 

The two-body contribution to the normalization integral is found from (21) 
or (21): 

= f»(0) = y^a / ■ r)!lj ' {±r lU l r • (23) 



(m 2 -x(l -x)M 2 Y 



In the limit of asymptotically small binding energy (a — > 0) we should substi- 
tute in (23) the function gM->2m(z), Eq. (13), and keep the leading term only 
(which now is ~ a = 1). In this case we obtain a trivial result: 

N 2 = 1. 

We see that the normalization integral is saturated by the two-body Fock sec- 
tor. The higher Fock sectors are generated by graphs with extra exchanges, 
each of them containing extra degree of coupling constant a. Their contri- 
butions are out of the leading a-term, which was only kept in the function 
9M->2m( z )- They are neglected in the limit a — > 0, that results in N 2 — 1. 

Now we calculate N 2 in next to the leading order, taking into account the cor- 
rection ~ a log a. For this aim, we substitute in the integral (23) the solution 
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(15). Calculating the integral (23) and keeping the terms ~ a log a, we obtain: 
In 1 

N 2 = l log-. (24) 

7r a 



Hence, for the contribution of the higher Fock sectors with n > 3: N n > 3 = 
1 — N 2 , we find: 

iV ri> 3 = ^logi. (25) 
7r a 



Below, in Section 6.2, we will compare Eq. (25) for iV n >3 = N$ + N4 + N$ + . . . 
with the three-body contribution N 3 . 



5.2.2 Extremely large binding energy (M — > 0) 

In this case we should substitute in (23) the function gM=o(z), Eq. (18), and 
put M = 0. This gives: 

} 9 
N 2 = F 2b (0) = 90 J x 3 (l - x) 3 dx = — w 64%. (26) 



The value of many-body contribution A^ n > 3 = 1 — N 2 at M = (B = 2m): 

00 

iy„> 3 = 1 - N 2 = 2 N n = 5/14 » 36% (27) 

n=3 

is the maximal value in the model. For M = the sum -^n over the 

many-body Fock sectors contains increasing degrees of huge coupling constant 
a = 2tt. However it still converges (to the value 5/14). 

5.3 Asymptotic at Q 2 — > 00 

5.3.1 Small binding energy (M — > 2m,) 

In the limit M — > 2m, form factor is obtained by substituting in (22) the 
solution (13) for <?M->2m(z)- It was calculated, by a different method, in [20]. 
The asymptotic formula reads: 



16a 4 m 4 

~Q 4 



F™ ymp {Q 2 ) 



2 - 



(28) 



14 



The limit M — > 2m corresponds, for massless exchange kernel, to a — > 0. 
As explained above, in the leading a order the contribution of higher Fock 
sectors is neglected. Therefore in this approximation the two-body contribu- 
tion coincides with the full one. They coincide not only in asymptotic, but 
everywhere. 



5.3.2 Extremely large binding energy (M — > 0) 



Substituting in (21) the expression (18) for gM=o(z) and making replacement 
of variable u(l — u) = v, we find at M = 0: 



1 1/4 

F 2b (Q 2 ) = 1080m 6 J x 3 (l - x) 3 dx J - 



v dv 



[m 2 + Q 2 vx 2 ) 3 y/l - 4v' 



(29) 



Form factor (29) is calculated explicitly: 



F 2b (Q 2 ) 



180m 4 



(Q 



2\3 



2m 



+3(Q 2 -20m 2 )Wl 



4m 2 . V4m 2 + Q 2 + ^ 



Q 2 



log 



2m 



At Q — > oo we find: 



F^ b symp (Q 2 ) 



270m 4 



M~ 2 



14 

y 



(30) 



(31) 



It is instructive to compare F2 b symp (Q 2 ) with full form factor and with the 
three-body one. This will be done below. 



6 Three-body contributions 



6. 1 Form factor 



Comparison of two-body contribution F 2 b{Q 2 ) with full form factor Ff u u(Q 2 ) 
(made below analytically for asymptotic domain and numerically for any Q 2 ) 
shows the total contribution to form factor of the higher Fock sectors with 
n = 3,4,5,.... In order to see the convergence of the Fock decomposition, 
we calculate in this section the first term of this sum, i.e., the three-body 
contributions. The diagrams are shown in figs. 8 (a) and (b). 
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Fig. 8. Three-body contributions to form factor. 

Figures 9 (a) and (b) also represent the three-body intermediate states. The 
graph Fig. 9 (a) incorporates self-energy, whereas the graph Fig. 9 (b) can 
be interpreted as a correction to the constituent form factor. However, in the 
ladder approximation, they do not contribute and therefore are omitted. There 
are no other three-body states. 




Fig. 9. Three-body states which do not contribute to form factor in the ladder 
model. 



The three-body form factor is calculated in Appendix C.2. The diagrams (a) 
and (b) in Fig. 8 give equal contributions. Therefore we take expression for 
diagram Fig. 8 (a) and multiply it by 2. The result has the form: 



F 3b (Q 2 ) 



am 2 r ip(Ri±, xi) ip(R'i±, x^) 9(x[ — x±) 
~2^J (<>) _ m 2 ) (4 a) - M 2 ) (A ~ x i) 
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d 2 R 1± d Xl d 2 ^ 1± d^ 

2xi(l-xi) 2xi(l -x^)' 1 ' 

The variables s^, are the energies squared in the intermediate states 2 
and 3 in Fig. 8 (a). They are expressed in terms of the integration variables 

— * — * — * 

R±, R'± and also Q± (see Appendix C.2): 



(a) _ R 2 1± + m 2 (Rfi^-Ri^ + x^Qif (R' 1± + x[Q ± ) 2 + m 2 

X\ X\ — X\ 1 — X\ 

(a) _ (jix ~ XxQ^f + m 2 (# 1± -R 1± + x 1 Q i _) 2 R\ L + m 2 

Xi — Xi 1 — 

We emphasize that in contrast to the two-body form factor (20), the variable 
R'i± in (32) is not expressed through R±± and Q±. Both R± and R'± are 
independent integration variables. 



6.2 Normalization integral 



Three-body contribution A 3 to full normalization integral is simply A 3 = 
F 36 (0), that with Eq. (32) gives: 



am 2 fil}(Ri±,Xi) if}(R'ij_,x' 1 ) 6(x' 1 -x 1 ) 

iV 3 — 



where 



/ 



2tt 5 J (s a - M 2 ) 2 (x; - x0 



d 2 /Z 1± dxi d 2 # 1± dx^ 
X 2xi(l-xi) 2x' 1 (l-x' 1 )' 1 j 



i? 2 ± + m 2 (R' 1± - R lA _f R'l ± + m 2 



+ ^ ^ + lliL^L. (36) 
Xi x x — Xi 1 — x x 



In the leading a log a order A 3 is calculated analytically in Appendix C.2. 2. 
The result is the following: 

A 3 = — log-, (37) 
7r a 



i.e., it coincides with A n > 3 , Eq. (25). It is natural, since the difference between 
A n > 3 and A3 is due to the higher Fock sectors A„> 4 , which contain extra 
degrees of a omitted in the present calculation. 
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The numerical calculations of N3 vs. M will be given in Section 8. At M = 
we get: N 3 = 0.257. This illustrates the convergence of the Fock decomposition 
at a = 2n: N 2 « 64%, N 3 « 26%, N 4 + N 5 + . . . « 10%. 



6.2.1 Asymptotic at Q 2 — > oo 

As discussed above, at asymptotically small binding energy only two-body 
contribution survives. Therefore we consider here the case of extremely large 
binding energy only. Since ip(R±,x) at large R± decreases like ~ 1/R± ( see 
Eq. (8)), the integral (32) well converges. Therefore we can take the limit 
Q 2 — > oo directly in the integrand. We put a = 2n, corresponding to M — 0. 
For the Q-dependent factor in (32) we find: 

1 _ (1 -x'){x' -xf 1 

Q ™oc ( S W -M 2 )(4 o) -M 2 ) ~ ^(1 - a;)*' 2 Q 1 ' 

Then the integral: 



1 x' 

1 rn 2 r , r (x' — x) 



Q 4 4tt 4 7 7 x 2 (l-a;)V 4 

v y 

xiI>m=o(R'±, x') d 2 R' ± ij M =o(R±,x) d 2 R ± 

with the wave function ■0m=o given by Eq. (B.5) is easily calculated and we 
get: 

FST* (Q 2 ) = (38) 



7 Full contribution 



7. 1 Form factor 



The contribution of all Fock components is incorporated by the form factor 
FfuiiiQ 2 ) in Eq. (9). Its calculation does not differ from the normalization 
integral explained in Section 3. The results is: 



1 1 1 1 

Ffuii(Q 2 ) = - ¥ ^ Jg M (2x-l)dx Jg M (2x'-l)dx' J u 2 (1 - uf 
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(xx'u(l - u)Q 2 + (6f - 5)m 2 + 2f(l - i)M 2 ) 
x - -4 du, (39) 

(xx'u(l - u)Q 2 + m 2 - £(1 - i)M 2 



where £ = xu+x f (l— u). At Q = form factor (39) turns into the normalization 
integral (10). 



7.2 Asymptotic at Q 2 — > 00 



As explained in Section 5.3.1, in the limit of asymptotically small binding en- 
ergy the higher Fock sector contribution disappears. Therefore the full form 
factor Ff u u(Q 2 ) coincides with the two-body one F 2b (Q 2 ), which has asymp- 
totic given by Eq. (28). 

We calculate asymptotic of Ff uU {Q 2 ) in the opposite case of extremely large 
binding energy. For this aim we substitute in (39) the function gM=o(z) from 
(18) and make the replacement v = u(l — u). This gives: 



1/4 



FfuiiiQ 2 ) = 540m 6 J V dV J x(l -x)dxj x'{l - x') dx' 

V 4 ~~ V 



fm 2 (5 - 3x - 3x') - Q 2 xx'v) 

x^ — '- (40) 

( m 2 + Q2 vxx iy ^ > 



The integrals over x and x' are calculated analytically: 



. 3 



1/4 

^ ,^9n 540m' ± f dv r ^9 / 9 ^9 n 

Fm(Q 2 ) = —gr- J ivr^ \- 2Q " (3m + e 



540m 4 /" df 

-(m 2 + Q 2 t;)(6m 2 + Q 2 V ) log ( 1 + ^] - 2m 2 Q 2 vU 2 ( 

\ m 2 J ' \ m z 

Now we take the limit Q 2 — > 00 and therefore keep in the integrand the leading 
term oc Q 4 - We obtain: 



3 h'i2(z) = Polylog(2, z) is a special function dilogarithm: 





Li 2 (z) = | log(l-t)j. 



The leading terms of Li2(— z) at 2 — > 00 are Li2(— z) « — i log 2 (z + 1) — 
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, V 4 
540m 4 r dv 



r dv 

J 



Q 4 J v^4^ 
o v 



2 - log 1 + 



QV 



It is also calculated analytically. Taking again the leading term at Q 2 — > oo, 
we find the asymptotic: 



F a„ (g2) 



270m 4 



(41) 



Comparing (41) with Eq. (31) for F 2 b, we see that at Q — > oo the leading terms 
27 q^ 4 log of Fj n // and F 2b , exactly coincide with each other. This means 
that asymptotic of Ff uU is dominated by F 2b and the many-body contributions 
decrease faster with increase of Q 2 than the two-body one. 

Comparing Eqs. (31), (38) and (41), we find the relation 

F fuu mP (Q 2 ) = F™ ymp (Q 2 ) + F™ ymp (Q 2 ) (42) 



valid for the terms ~ 1/Q 4 - We see that the term ~ i log ^ in Fj^ np comes 
from F™ ymp and the difference between F^™ p and F^ ymp Q f the order of 
~ ^4 results from F^ ymp (Q 2 ). The corrections to Eq. (42) result from the 
higher contributions F^, F 5b , etc. Since they are absent in (42) in the order 
^4, they decrease more rapidly than This shows a hierarchy of asymptotic 
form factors: F 2b > F 3b > F ib > . . .. This also indicates a good convergence 
of contributions of the Fock components with increasing particle numbers, a 
least, in asymptotic. In Section 8 we will calculate numerically the different 
contributions to the normalization integral (i.e., take the opposite, Q 2 = 
limit) and also find good convergence. 



7.3 Arbitrary binding energy 



The momentum transfer Q 2 enters the equation (39) for form factors Ff u u in 
the combination xx'Q 2 , and in the equation (21) for F 2b as x 2 Q 2 . Therefore, 
the form factors at Q 2 — > oo are determined by the integration domain near 
x, x' — > and, hence, by behavior of the function gu{z = 2x — 1) at binding 
point z — — 1. As shown in Section 4, at z — > — 1 the function <?m(z) linearly 
tends to zero: <?m(z ~~ * — 1) ~ (z + 1), independently of the value of M. 
Therefore, the asymptotic F fuU (Q 2 ) ~ F 2b (Q 2 ) ~ ^ log(£a ) is expected for 
any M. This is natural, since this asymptotic is valid in both limiting cases 
M = and M w 2m. 
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To check that, we represent Qm{z) as: Qm{z) = {z 2 — l)f(z), where the factor 
(z 2 — 1) provides zeroes at z — ±1 and f(z) is an even function of z which can 
be decomposed near z = as: f(z) = c + c 2 z 2 + c 4 z 4 + . . .. When we take the 
zero degree in this decomposition, i.e. the term c only, we come back to the 
case gu{z) oc (1— z 2 ), considered in the case M — 0, and therefore for arbitrary 
binding energy we obtain the asymptotic F fuU (Q 2 ) ~ F 2b (Q 2 ) ~ ^log^) 

found for M — 0. Other terms of decomposition result in the factor z 2m z' 2n = 
(2x — l) 2n (2x' — l) 2n which does not change the leading term. The coefficient at 
the leading term and next to the leading terms depend on particular function 

/(*)■ 



8 Numerical calculations 

The contributions N 2 and N 3 to the normalization integral from the two- 
and three-body sectors are given by Eqs. (23) and (35) correspondingly. The 
difference: 

iV n > 4 = 1 - N 2 - N 3 

is contribution to full normalization integral of higher Fock sectors (i.e., prob- 
ability of the sectors with n > 4). 

These probabilities, i.e., N 2 , N 3 and iV n >4 as function of M are shown in 
Fig. 10. Mass M (and momentum transfer Q below) are given in units of m, 
i.e., we put m — 1. The maximal values of 7V 3 and N n > 4 are achieved at M — 0. 
In this case N 2 is calculated analytically by Eq. (26): N 2 = 9/14 w 0.643. A^ 3 
is found numerically: A^ 3 m 0.257. Hence, the total contribution of the Fock 
sectors with n > 4: N n > 4 (M = 0) = 1 - N 2 - N 3 = N 4 + N 5 + . . . w 0.100 
is 10% only. Like in the case of form factor, this shows quick convergence 
relative to increase of the particle number in the Fock sectors. In the interval 
< M < 1.8 these contributions are rather smooth functions of M and then, 
when M tends to 2, N 2 tends to 1 and both A^ 3 and N n > 4 tend to zero very 
quickly. 

However, N n > 4 decreases faster than jV 3 . This is seen from Fig. 11, where the 
ratio N "r* vs. M is shown. This ratio decreases at M — > 2 and N 3 dominates 
over N n > 4 . This results from the fact that the higher Fock sectors iV n >4 contain 
extra degree of a. This is also reflected in the coincidence of analytical formulas 
(25) for iV„> 3 with (37) for N 3 . 

In Fig. 12 the form factor Ff u u(Q 2 ) is shown vs. Q 2 for the M values M — 0, 
1.5, 1.9, 1.99, 1.9999. It is calculated by Eq. (39). We emphasize very strong 
dependence on M. For example, at Q 2 = 1 the form factor for M = 1.9999 is 
many orders (~ 10 6 times) smaller than the form factor for M — 0. 
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Fig. 10. Contribution to the full normalization integral of the states with the con- 
stituent numbers n = 2, n = 3 and n > 4. 




0.2- 



0.1 - 



0.0 



0.0 



0.5 



1.0 



M 



1.5 



2.0 



Fig. 11. Ratio ^± vs . m. 



In order to show, how the two-body form factor achieves its asymptotic value, 
in Fig. 13 we present the ratios F™ ymp (Q 2 ) / F 2b (Q 2 ) for M = and M = 
1.9999. For M — > 2m the asymptotic form factor F^ 6 S2/mp ((5 2 ) is given by 
Eq. (28), and for M = - by Eq. (31). We see that for small binding energy, 
M = 1.9999, the asymptotic is achieved much earlier (at Q 2 ~ 0.1m 2 ) than 
for large binding energy M = (Q 2 ~ 1000m 2 ). 
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Fig. 12. Form factor Ff u u(Q 2 ) vs. Q 2 for a few values of M. 




Fig. 13. Ratio of asymptotic of form factors F% b symp (Q 2 ) / F^Q 2 ) ■ 

Similar situation takes place for the ratio of full form factors F^ l mp (Q 2 ) /Ff u u(Q 2 ). 
As shown in Section 5.3.2, at large values of Q 2 the many-body contribution 
decreases more rapidly than the two-body one. Therefore at Q 2 — > oo the 
form factor Ff uU {Q 2 ) coincides with F 2 b(Q 2 )- The ratio of the form factors 
F 2 b(Q 2 )/F fu u(Q 2 ) vs. Q 2 for M = 1.9999 and M = is shown in Fig. 14. 
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Fig. 14. Ratio of form factors F 2b (Q 2 ) / 'F fuU (Q 2 ). 



For very small binding energy B = 0.0001 (M = 1.9999), and, correspond- 
ingly, small a, the many body contribution is practically absent. Therefore 
F 2b (Q 2 ) I 'Ff u u(Q 2 ) a* 1 for any Q 2 . For very large binding energy B = 2 
(M = 0) the ratio F 2b (Q 2 ) / 'FfuiiiQ 2 ) tends to 1, but it becomes more or less 
close to 1 at huge values Q 2 ~ 10 20 m 2 , when log(Q 2 /m 2 ) dominates. Both form 
factors F 2b {Q 2 ) and Ff u u(Q 2 ) contain \og(Q 2 /m 2 ) with the same coefficient 
but they differ by the terms without \og(Q 2 /m 2 ) (the terms —14/3 = — 4| 
in F 2b (Q 2 ) in (31) and -4 in F fuU (Q 2 ) in (41)). The difference of form fac- 
tors is — |. Therefore the form factors F 2b {Q 2 ) and Ff u u(Q 2 ) become to be 
close to each other very far, when log(Q 2 /m 2 ) ^> 2/3. The relative difference 
(FfuiiiQ 2 ) - F 2b (Q 2 ))/F full (Q 2 ) decreases as ~ l/log(Q 2 /m 2 ) only. 

However, in a limited domain of Q 2 the form factors Ff u u(Q 2 ) and F 2b {Q 2 ) 
differ approximately by a factor. It is illustrated by Fig. 15 for M = 1.5. Solid 
line in this figure represents the form factor Ff u u(Q 2 ). Dash line is F 2b (Q 2 ). Its 
value at origin is: N 2 = F 2b (0) = 0.669 (compare with N 2 = F 2b (0) = 9/14 = 
0.643 for M — 0). The dotted curve shows the two-body form factor F 2b (Q 2 ) = 
F 2b (Q 2 ) I F 2b (0) , normalized to 1 at Q 2 = 0. It is almost indistinguishable 
from Ff u u(Q 2 ). If this coincidence takes place in a more realistic model, then, 
comparing the experimental data with the two-body form factor, normalized 
to 1 (as it is usually done), we would conclude that the later dominates and we 
would not notice 33% many-body contribution (N n > 3 = 1 — 0.669 = 0.331). 
This would be wrong conclusion about the structure of a system. Therefore 
the results shown in Fig. 15 could be instructive. 
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Fig. 15. Form factors for M = 1.5. Solid line - F fuU (Q 2 ). Dash line - F 2b (Q 2 ). 
Dotted line - F 2b (Q 2 ) = F 2b (Q 2 )/F 2b (0). 

9 Summary and discussion 



The state vector in field theory does not correspond to a definite number of 
particles. It is represented as a superposition of Fock sectors with different 
particle numbers n. In Wick-Cutkosky model this sum starts with n = 2 
and goes until infinity. We studied the contribution to the electromagnetic 
form factor (and, in particular, to the normalization integral) of the Fock 
sectors with n = 2, n = 3 and sum of all contributions with 2 < n < oo. 
For small and large binding energy (and coupling constant a) the results are 
obtained analytically (except for the three-body contribution at M = 0), and 
for arbitrary binding energy - numerically. 

In the limit a — > the two-body sector survives only. Corresponding two-body 
Fock component is simply the non-relativistic ground state wave function in 
Coulomb potential. 

In the next order ~ a log a the higher Fock sectors contribute, however the 
two-body sector still dominates. The correction ~ a log a comes from the 
three-body contribution, the correction ~ a 2 results from many-body sectors 
with n > 4. Since at small a the higher Fock sectors are suppressed by extra 
degrees of a, the Fock decomposition converges very fast. Not only the nor- 
malization integral, but the form factor Ff u u(Q 2 ), at any Q 2 , is also dominated 
by the two-body contribution. 

In the leading order a log a the contribution of many-body sectors to the nor- 
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malization integral N n > 3 , extracted from the BS function, Eq. (25), coincides 
with the three-body contribution A3, Eq. (37), found by direct calculation 
in LFD of the three-body intermediate state. This coincidence illustrates the 
main idea of our work: the "two-body" BS function indeed includes many- 
body Fock sectors. Therefore both the two-body contribution N 2 and sum 
of many-body contributions A n > 3 can be found from a given "two-body" BS 
function. 

In the opposite case of large coupling constant, up to its maximal value a = 
2n (corresponding to M — 0) the Fock decomposition still well converges. 
Contributions to the normalization integral of a few Fock sectors is presented 
in the Table 1. In spite of huge value of the coupling constant, contributions 
of the Fock sectors quickly decrease with increase of number of particles. Two 
first sectors with n = 2 and n = 3 determine 90% of the normalization integral. 



N 2 


N 3 


iV„> 4 


N 2 + N 3 + N n > 4 


0.643 


0.257 


0.100 


1 



Table 1 

Contributions of the Fock sectors with the particle numbers n = 2, n = 3 and n > 4 
(N n >4 = YHn=A N n ) to the full normalization integral N = J2™=2 N n = 1 of the 
state vector for M = (a = 2ir). 

This can be explained as follows. Any extra Fock sector adds to the ampli- 
tude extra factor ~ a/(s — M 2 ). There is a competition between increasing 
coupling constant in numerator and increasing energy in denominator. When 
a increases, the binding energy increases too. Then the bound state mass M 
decreases. This results in increase of the energy denominators (s — M 2 ) so 
that effective contribution of the factor ~ a/(s — M 2 ) remains smaller than 
1. This is a probable mechanism of suppression of the higher Fock sectors in 
the Wick-Cutkosky model. 

The behavior of form factors vs. Q 2 is also rather instructive. The decrease of 
full form factor Ff u u(Q 2 ) (Fig. 12) with increase of Q 2 strongly depends on 
the bound state mass M. At smaller values of M the form factor begins to 
decrease at larger values of Q 2 . At M = the form factor is almost constant 
up to Q 2 10m 2 and decreases at larger Q 2 . 

Asymptotic of form factor for the larger binding energy B is achieved at larger 
values of Q 2 (Fig. 13). For mass M = 1.9999 form factor F 2 b{Q 2 ) obtains its 
asymptotic form at Q 2 > m 2 , whereas for M = - at Q 2 > 1000m 2 . 

The contribution of Fock sectors to form factor decreases with increase of 
Q 2 more faster for higher Fock sectors (for larger values of n). Therefore, 
in asymptotic only the two-body contribution survives (Fig. 14). For small 
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binding energy (M = 1.9999) the two-body component dominates for any Q 2 . 
With increase of binding energy the many-body components become more 
important. However, they are suppressed at larger values of Q 2 . The value 
of Q 2 when many-body contributions can be neglected is larger for larger 
binding energy. For M = the many-body contributions become smaller than 
5% (and the two-body form factor constitutes 95% of full form factor) at very 



The leading asymptotic terms ~ ^ log ^ in Ff u u and F 2 b are the same (with 
the same coefficient). Therefore, at Q 2 — > oo: Ff uU = F 2 b- The form factor F 3b 
does not contain log but it contains next to the leading term ~ With 
this correction Ff u u at Q 2 — > oo is determined by sum of two form factors, 
Eq. (42): Ff u u = F 2 b + Fsb, where F^ decreases faster than Fib- Contributions 
F n >4 of the Fock sectors with n > 4 is the next correction which decrease 
faster than F 3b - 

Though form factor Fib substantially differs from Ff u u, in a limited but enough 
large domain of Q 2 their ratio is a constant. Being normalized at Q 2 = to 
1) F 2 b{Q 2 ) becomes to be very close to Ff u u(Q 2 ) everywhere (see Fig. 15). 
Approximate description of true form factor of a system by the two-body 
one (with artificially imposed condition F 2b (0) = 1) does not mean that this 
system is the two-body one. 

Good convergence of the Fock decomposition in Wick-Cutkosky model in- 
dicates that the Fock decomposition may be efficient in more realistic field 
theories, which, however, should be studied separately. 
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A Correction ~ a log a to the Wick-Cutkosky solution 

Eq. (11) with the kernel (12) is equivalent to the following differential equation: 



large Q 2 > 10 10 m 2 . 



9m(z) + 



a 



m 2 g M (z) 



= 



(A.l) 



7T (1 - z 2 ) m 



2 - \(l-z 2 )M 2 
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with the boundary conditions ^(±1) = 0. In the limit a — > its solution (13) 
does not depend on a (except for a normalization factor). To find a correction 
to the Wick-Cutkosky solution (13), we solve Eq. (A.l), using the method of 
the paper [14]. Namely, the Eq. (A.l) is identically rewritten as: 

gl(z) + HVo + V 1 )g M (z) = 0, (A.2) 



where 



1 1 2 4m 



rj 



M 2 



1, V = S(z), Vi = - 



IT 



l+rf 



1 - z 2 )(r] 2 + z 2 ) 



8{z). 



Since the lowest order solution (13) g^'{z) = gM->2m{z) ~ (1 — \z\) is obtained 
from (A.2) at Vi = 0, i.e., it satisfies the equation g" (0 \z) + \ (0) V g^(z) = 
with A(°) = 2, Vi can be considered as a perturbation. Eq. (A.2) becomes: 



g 



//(0) , „//(!) 



+ g" w + ... + (A<°> + A« + • • .)(Vo + + «,<U + ...) = 0, 



which in the first order gives [14]: 



f\gWV ig W dz 



and the simple equation for g^: 

'\ i0) V g^ + A<°)Vi<7< ) + A (1) V~ # (0) 



(A.3) 



From the relation A = A^ + \^ = - with A^ = 2 we find the binding energy 

(16) and then the solution g^(z) of Eq. (A.3). It contains the terms ~ a and 
~ a log a. We have found that the next order correction g^ 2 \z) (relative to the 
order of the kernel Vi) contains the linear a contributions (i.e., the same order 
of a as appears in g ( - 1 \z)). But it does not contain the terms ~ a logo;. The 
binding energy calculated in [14] by another method (in LFD) coincides with 
(16) up to the ~ a log a: terms too. Therefore we keep the leading ~ a log a 
term only. In this way we find eg. (15) for gui 2 )- 



B Explicit BS solutions and the LF wave functions 



Substituting (13) into (7), we find the BS function for asymptotically small 
binding energy explicitly: 
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-4i v / 27TQ; 5 / 2 m 3 
= ^ (B.l) 

x 7 r , (B.2) 

(m 2 - (\p + A;) 2 - iO) (m 2 - (|p - A;) 2 - io) 



where M = 2m - 



Substituting (18) into (7), we find the BS function for extremely large binding 
energy: 



„ , 6i^nm 3 (k 2 - m 2 ) 

<f> M =o(k,p)- 



(p-k) 2 (m 2 — k 2 — p-k — iO) (m 2 — k 2 + p-k — iO) 

3i\/307rm 3 ^ m 2 — k 2 — p-k — iO 
(p-k) 3 ^ m 2 — k 2 + p-k — iO 



By means of Eq. (8) we find the LF wave functions for these two limiting 

cases: 

, T x 8V^a 5/2 m 3 x(l - x) \ x, if < x < \ 

^M->2 m (k±, x) = — - ^ x 2 (B.4) 

(kl + m 2 -x(l -x)M 2 ) [ 1 -x, if \ < x < 1 

t 12V3dnm 3 x 2 (l-x) 2 
■^M=o{k ± ,x) = — -2 . (B.5) 

(k\ + m 2 ) 



C Calculating form factors 



C. 1 Two-body form factor 



To calculate the LF amplitudes, we use the Weinberg rules [11]. To find the 
amplitude —At, one should put in correspondence: to every vertex - the factor 
g (in the theory with interaction H mt = —gip 2 (x)x(x)), to every intermediate 
state - the factor: So _ 2 s+i0 , where s is the energy squared in the intermediate 
state, Eq. (3), and so is the initial (=final) state energy. In our case (bound 
state): s 2 , = M 2 . To every internal line one should put in correspondence the 
factor: One should take into account the conservation laws for R i± and 
Xi in any vertex and integrate over all independent variables with the measure: 
ro \z ■ The variables R i± introduced by Eq. (2) in Section 2. 
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The LF graph determining the EM vertex of two-body state is shown in Fig. 7. 
The energies squared corresponding to the intermediate states 1 and 2 in Fig. 7 
read: 



B - Or MlrV- R2 l± +m2 , R 2± + ^ _ gl±+^ 

Si — [Ki + K 2 ) — 1 — 7- r, 



■ 2 



- = <*' + k ^ = §J0-y (ai » 

We take into account: R[ = k± — x±p' — R\ — x±Q, where R± = k\ — X\p and 
Q = p' — p. Hence: R' 1± = Ri± — XxQ±. Using the rules [11], we get: 



(p + p') P F 2b = J(k 2 + k' 2 ) p r 



M 2 — s' M 2 -s 



y 9(x 1 )9(x 2 ) 6(x 2 ) d 2 R 1± d Xl 

X 2xi 2x 2 2x 2 (2vr) 3 ' [ ' } 

We used the standard condition ou-Q = 0. Therefore u-k 2 = u-k' 2 and, hence, 
x' 2 — x 2 — 1 — x\. T is the vertex function related to the wave function as: 

* = lhn- (a3) 



Multiplying both sides of the equality (C.2) by u p , we obtain the formula (20). 



C.2 Three-body form factor and normalization integral 



C. 2. 1 Form factor 

We find form factor F^f,, calculating directly the amplitudes of diagrams Fig. 8. 
This is of course equivalent to perturbative calculation of the three-body com- 
ponent and then expressing form factor through it by a standard formula. 
Using the Weinberg rules [11], similarly to two-body contribution, Eq. (C.2), 
we obtain for the three-body graph Fig. 8 (a): 



(p + p) P F$=g 2 J(h + k' 2 ) p 

XT' 2 2 2 2 

(M 2 -4 a) ) (M 2 -4 a) ) (M 2 -4 a) ) (M 2 -4 a) ) 
9(x' 1 ) 9(1 - x\) 9(1 - x\) 9(x' l - Xl ) 

X 2x[ 2(1 - xi) 2(1 - xi) 2(xi - xi) 
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9( Xl ) 0(1 -m) tfR^dx[ d 2 R 1± d Xl 
X 2xi 2(1 -xi) (2tt) 3 (2vr) 3 ' 1 ' ' 

Here we used that x 3 = x\ — x±, x& = x' 2 = 1 — x[, x-i = 1 — x±. The 
denominators (M 2 — s^) and (M 2 — s^) are absorbed in the relation (C.3) 
between T and ip. The values of s 2 a \ s 3^ are defined in (33), (34). Multiplying 
both sides of Eq. (C.4) by uj p , we obtain the contribution of graph Fig. 8 (a), 
which is half of Eq. (32). 

Expression (33) for is found as follows. Since in the intermediate state 
2 in the graph 8 (a) there are three particles (with the momenta ki, k 3 , k 4 ), 
according to Eq. (3) we have: 



4 a) = (h + h + k A f = ^ + m2 + ^ + Rl +m2 



Xi x 3 x 4 

where R\±, R-3±, R 4 ± are the perpendicular to uj components of the four- vectors 

Ri — k\ — xip, R 3 = k 3 — x 3 p, R4 — k 4 — x 4 p. 

Using the momenta conservation in the vertices, we transform R 3 as: 

R3 = k 3 — x 3 p = k[ — ki — (x[ — xi)p (mod uj) = R[ — R 1 + x[Q, 

where (mod uj) means that we omit the terms, proportional to uj since they 
do not contribute to the _L-components. The variable R[ is defined as R[ = 

— * — * — * — * 

k'i—x'ip'. Hence we get: R 3 ± = R'i± — Ri±+x' 1 Q±, as appears in (33). Similarly 
we find: R 4 ± = —R'i± — x[Q± and reproduce in this way Eq. (33) for s^. 

To find s 3 a \ we start with: 

~J2 ~ 2 _ o 



4 a) = ( h + k 3 + k> ) 2 = + ^ + ^li±^ , 



Xi X 3 X 

where all i?'s are defined relative to p'\ 

Ri = ki- x x p\ R 3 = k 3 - x 3 p', R[ = k\ - x\p' 

and then, again using the conservation laws, we find R±± = R±± — XiQ±, 
R 3 ± = R\± — R±± + XiQ± and derive Eq. (34) for s 3 a \ 

Similarly we obtain the contribution of the diagram Fig. 8 (b): 
l(6) 167ram 2 f 4>(Ri±,xi) ij(R'i±,x[] 



ft. 



3b (2vr)6 /(4 6 )_M 2 )(4 b) -M 2 ) 
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e{x l -x , 1 ) d^gixd^i d 2 # 1± dxj 

X (xi-xi) 2x 1 (l-x 1 ) 2^(1 -xi)' ^ ' ^ 



where 



(&) _ (i?' 1± +xjQ ± ) 2 + m 2 (E 1± - R' 1± - xjQ ± ) 2 R 2 ± +m 2 
s 2 — / 1 / 1 ; 

X\ X\ — X\ 1 — X\ 

(b)_R'l±+m 2 (R 1± - R' 1± - Xl Q ± ) 2 (R 1± - Xl Q ± ) 2 + m 2 
s 3 — 1 1 ■{^■') 

X 1 X\ — X\ 1 — X\ 



By the replacement of variables R\± — > R'i±,x\ — > x[, R\± — > i?i_i_,xi 
xi we show that F^j = F$ ■ Therefore we finally get for full three-bc 
contribution: 

F 3b (Q 2 ) = F^\Q 2 ) + F£\Q 2 ) = 2F^(Q 2 ) 
and derive in this way Eq. (32) for F 3b (Q 2 ). 



C.2.2 Normalization integral 

Three-body contribution to the normalization integral N 3 = F 3fe (0) is given 
by Eq. (35). It is instructive to represent it through the interaction kernel. For 
this aim, we rewrite (35) as 



N 3 = Ft\0) + F^>(0) 

,2 



with 



am 



J ^(R'i±,x[ 



Air 5 
x^(-Ri±,xi) 



9(x[ - Xl ) 



+ 



0(xi-a;i) 



d 2 R 1± dxi d 2 R' 1± dx\ 



2xi(l - xi) 2xi(l - xi) 



(Ci 



«6 = 



i?V + m 2 



| (^i±-^i±) 2 | 



i? 2 ± + m 2 



Xi 



Xi 



Xi 



1 — Xi 



(C.9) 



The ladder kernel is shown graphically in Fig. 16. Its calculation by the Wein- 
berg rules [11] gives: 

Area 6>(xj - xi) Ana 6(x 1 - xj) 

(xi-xO^-M 2 ) (xi-xi)( S6 -M 2 )' ^' Uj 
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Fig. 16. Ladder kernel. 



with s a ,Sf, given by Eqs. (36), (C.9). The kernel (CIO) (which enters the 
Schrodinger-type equation) is related to the amplitude /C of graph Fig. 16 as 
V = -/C/(4m 2 ). 



Let us calculate the derivative dV/dM 



2. 



dV 
dM 2 



-Ana 



6{x[ - Xl ) 



+ 



(xi - x l )(s a - M 2 ) 2 ( Xl - x'OK - M 2 ) 2 



(C.11) 



Comparing (C.ll) with the integrand of (C.8), we represent N 3 , Eq. (C.8), in 
the form 



d 2 R 1± dxi d 2 i?' 1± d^ 

x ; r — — tt- (C.12) 

2^(1 -xi) 2x' 1 (l-x' 1 ) v ; 



Then we make replacement of variables R±±,Xi — > k, 9: 

R 1± = ksm9, Xl = l(l-^l), (C.13) 

2 V £ k J 



where e k = Vk 2 + m 2 , and similarly for R' 1±: x' 1 — > k',6'. The integration 
volume is transformed as: 

d 2 i?i_L dxi Ri± dR 1± dxid(/> k 2 dk sin OdOdcj) d 3 k 



2a;i(l-xi) 2xi(l-xi) e k e k 
In these variables N 3 obtains the form: 

Am 2 f ,.-> . dV(k',k,n,M 2 ) . d 3 A; d 3 ^ . 

N * = -(2^J^ k ' n) OM 2 ^ ~ ' (C ' 14) 



where we introduced the vector n: k-n = kcos9, k'-n = k' cos 9'. The equation 
(C.14) coincides with the many-body contribution (for a general kernal V, not 
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only for the ladder one), given by Eq. (3.50) in [10]. It can be obtained by 
more general method [10], without precising the kernel V. 

For small a the value N3 can be calculated explicitly. We rewrite N 3 , Eq. (35), 
as: 



am 2 rip(Ri±,x 1 )il>(R , 1 ±,x' 1 )(x' 1 -x 1 ) 

VV3 



/ 



2tt5 J [(s a — M 2 )(x' 1 — Xi)] 2 

, d 2 R 1± dx 1 d 2 R' 1± dx[ 

x9(x \ - xi) - — r - . — r- (C.15) 

v 1 ; 2x 1 (l-x 1 ) 2x' 1 (l-x[) v 1 

with s a represented more explicitly: 

R 2 1± +m 2 R' 2 1± -2R 1± R' 1± cos (j) + R 2 ± R\ ± + m 2 

Sa = + } + —, — ■ (C.16) 

X 1 X\ — X\ 1 — x\ 



(j) is angle between Ri± and R'i±. We make replacement of variables by Eqs. 
(C.13). Since N3, Eq. (C.15), is proportional to a, for the wave function we 
can take Eq. (B.4) for asymptotically small a, where higher degrees of a are 
neglected. In the variables (C.13) it obtains the form (cf. Eq. (3.64) from [10]): 

<KM) = 8V 7* 5/2 fc , wr , (C17) 



where k = \JmB = ma/2 (B is the binding energy). 

Then we make another replacement of variables k, k' — > p,p': 

k = up, k' = up, (C.18) 

where p,p' are dimensionless. In the limit k — > 0, we decompose the factor 
(s a — M 2 )(x[ — xi) in the denominator of (C.15) in series of k, keeping the 
leading term: 



(s a - M 2 )(x[ -xi)k,k 2 \p 2 - 2pp cos 9 cos 9' - 2pp sin 9 sin 9' cos + p 



/2 



K 



H (p +p + 2)(pcos9 -p'cos9') 



m 



(C.19) 



Introducing vector q = p' — p, we represent (C.19) as: 



(s a - M 2 )(x[ - xi) w n 2 (q 2 --(p 2 + p' 2 + 2)q-n). (C.20) 
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The difference (x\ — x\) in the numerator of (C.15) is transformed as: 

x\ — x\ ~ — (pcos^ — r> cos 9') = q-n. 

1 2m yF F ' 2m H 



Wave function (C.17) obtains the form: 
8^m~K 5 / 2 

«p) (c ' 21) 



It does not depend on 0, this simplifies the calculation. The integration volume 
reads: 

d 3 k k, 3 d 3 p 

r*~> 

Hence, TV 3 is transformed as: 

32nm f (n-q) 6(-n-q) p 2 dpdo^p' 2 dp'dop 

TV 3 — - — / 2 • (L<.22j 

71 J (mq 2 - k( P 2 + p' 2 + 2) n-<f) (p 2 + l) 2 (p' 2 + l) 2 



One can substitute here denominator in the form of (C.19) and integrate 
over three angles 6, 9' and <fi. However, the same result can be found in a 
more simple way. Namely, since TV3 in (C.22) does not depend on n, but 
the integrand depends on n, one can at first average integrand over the re- 
directions. Therefore we replace the factor: 

(n-q) 6 (—n-q) 

hi = 



(mq 2 — n(p 2 + p 12 + 2) n-qj 

by 



ho 



(n-q) 6 (—n-q) do?, 



(mq 2 - n(p 2 + p' 2 + 2) n-q) AlT 

The angle integration here is one-dimensional and h 2 is easily calculated. After 
that the angle integrations in (C.22) are also carried out easily. Then the 
integrals over p and p' are calculated approximately, using that a C 1 implies 
k < m. In this way we obtain the analytical formula (37) for TV 3 at a <C 1. 
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